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Abstract — In this work, the governing differential 
equations of elastic column buckling represented by a 
system of three coupled differential equations in the three 
unknown displacement functions, v(x), w(x) and 6(x) are 
solved using the method of Fourier series. The column 
was pinned at both ends x = 0, x = l. The unknown 
displacements were assumed to be a Fourier sine series of 
infinite terms, which was found to satisfy apriori the 
pinned conditions at the ends and substituted into the 
governing equations. The governing equations were found 
to reduce to a system of algebraic eigenvalue - 
eigenvector problem. The buckling equation was found to 
be a cubic polynomial for the general asymmetric 
sectioned column. The buckling modes were found as 
flexural torsional buckling modes. For columns with 
monosymmetric sections, it was found that the buckling 
mode could be flexural or flexural - torsional depending 
on the root of the cubic polynomial buckling equation 
which is the smallest. For columns with bisymmetric 
sections, it was found that the buckling modes are 
uncoupled and bisymmetric columns could fail by pure 
flexural buckling about the axes of symmetry or pore 
torsional buckling. The findings are in excellent 
agreement with Timoshenko ’s solutions. 

Keywords — Monosymmetric columns, bisymmetric 
columns, flexural torsional buckling mode, algebraic- 
eigen vector problem. 
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I. INTRODUCTION/LITERATURE REVIEW 

Euler’s [1] work on column flexural buckling presented 
the first analytical method of determining the buckling 
strengths of slender columns. St Venant [2] later worked 
on uniform torsion and presented the first reliable work 
on the twisting response of structures to torsion. Flexural 
torsional buckling was studied by Michell [3] and also by 
Prandtl [4] who considered the lateral buckling of beams 
of narrow rectangular cross sections. 

Prandtl and Michell’s work were extended by 
Timoshenko [5] to include the effects of warping torsion 
in T section beams. The works of Wagner [6], Vlasov[7], 
Timoshenko [8] and Timoshenko [9] led to the 
development of a general theory of flexural torsional 
buckling. Other researchers who have studied flexural 
torsional buckling of structures include Nwakali[10], 
Timoshenko and Gere [11], Alsayed[12],Trahair[13], 
Zlatko[14] and Al-Sheikh [15]. 

n. THEORETICAL FRAMEWORK 

The generalized elastic column buckling problem 
formulated in terms of the displaced configuration is 
represented by the following set of three linear 
differential equations obtained when the non linear terms 
are neglected [ 16] [ 1 3] . 
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where Mj z , Mj y — moments caused by transverse loads only. 

Ij 7 = Iyy + I— + (e 2 + e 2 )A = polar moment of inertia about the shear centre. 
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/ = St Venant torsional stiffness of the section 

E = Young’s modulus of elasticity 

G = shear modulus 

Iyy, I zz = moments of inertia 

C w = warping constant 

Iyz = product of inertia 

e y , e z = coordinates of the shear center 

V y , V z = shear forces 

q y , q z = transverse loads 

N x = axial load 

v(x),w(x) are transverse displacements 
9v = twist rotational displacement 

x = longitudinal axial coordinate, yz is the plane of the cross section 
A = area of the cross section of the column. 

The governing equilibrium equations for the generalized column buckling problem is represented by the system of three 
simultaneous differential equations in the three unknown displacements v(x), w(x) and 9(x). For columns with prismatic 
cross-sections, the elasticity properties (GJ and E) as well as the inertial and geometrical properties are constant; but the load 
coefficients (My, My, V y , V z , W w , V w ) are variables depending on the axial longitudinal coordinate (x). The system of 
governing differential equations thus have variable coefficients, rendering them difficult to solve mathematically. 

However, simplifications of the system of governing differential equations can be obtained by using special characteristics of 
different types of column problems. A simplification of the governing equilibrium equations can be obtained if the yz 
coordinates are principal coordinates. Then I yz = 9, and the governing differential equations become [ 16] [ 1 3] 



(2.4) 


(2.5) 


www.ijaers.com 


Page I 293 


International Journal of Advanced Engineering Research and Science (IJAERS) 
https://dx.doi.Org/10.22161/iiaers.4.3.45 


[Vol-4, Issue-3, Mar- 2017] 
ISSN: 2349-6495/P) 1 2456-1908(0) 


-(M lz -e y N x ) 


d 2 w 


dx 


+ v z *-v*L- 

2 ' dx dx 


r H V ^ 

r V +C V + w w 

'-'z v y ^ ^y v z ^ s, 

\ y 


de , A 

— = t(x) 
dx 


(2.6) 


If the member is free of transverse loads q y , q : become zero; and the transverse moments are constant, then the shear forces 
V y , V- are also zero. For member (columns) subjected only to axial compressive load N x acting through the centroid of the 
section, the moments due to the transverse loads vanish, the applied torque vanishes; and if the load is applied such that the 
bimoment vanishes, the system of differential equations become [15][16] 
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The system of differential Equations (2.7), (2.8) and (2.9) represent the governing equilibrium equations for a particular case 
of the generalized elastic column buckling problem presented by Equations (2.1), (2.2) and (2.3) when I yz — 0, 

q y = q- = 0, V y = (E = 0, M Iy = M Iz = 0, W w = 0 . 


III. APPLICATION OF THE FOURIER SERIES 

METHOD 

We seek to apply the Fourier series method to solve the 
system of differential Equations (2.7), (2.8) and (2.9) 
representing the governing-equilibrium equations for the 
elastic column buckling problem for the case when the 
ends at x = 0, x = /are on pinned supports. 

For pinned ends at x = 0, x = l, the relevant boundary 
conditions are 

v(x = 0) = 0 v"(x = 0) = 0 

w(x = 0) = 0 w"(x = 0) = 0 

0(x = 0) = 0 0"(x = 0) = 0 (3.1) 

v(x = /) = 0 v"(x = /) = 0 

w(x = /) = w"(x = /) = 0 

0(x = l) = 0"(x = /) = 0 
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Suitable Fourier series representations of the three 
unknown displacement functions that satisfy apriori the 
above boundary conditions are given by 
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where v m , w m and 0,„ are infinite number of unknown 
coefficients of the Fourier sine series representations of 
the unknown displacement functions v(x), w(x) and 9(x) 
that we seek to determine. If Equation (3.2) represent 
solutions of the governing Equations (2.7), (2.8) and 
(2.9), then 
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The stability problem represented by the system of differential equations is now reduced to an algebraic eigenvector 
eigenvalue problem given by the homogeneous Equation (3.12). For non-trivial solutions, the characteristic buckling 
equation is the determinatal Equation (3.13) 
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This is the characteristic equation for determining the 
buckling load of a column with an asymmetric cross- 
section. The buckling equation is a third order polynomial 
in N x . Thus, it can be solved using the methods for 
solving polynomials to obtain the three roots P cr ^ , P cr ^ 

and P cr ^ . The smallest of the three critical buckling loads 

will govern the buckling behaviour of the column. 

However, if the zz axis is the axis of symmetry of the 
cross-section, e y = 0, and the characteristic buckling 
equation simplifies to become 
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Hence, in this case (monosymmetric cross-section), the 
column can buckle in two possible buckling modes 
namely pure flexural buckling in the yy - axis, and 
flexural torsional buckling, otherwise. If Pyy is the 
smallest of the three roots, of Equation (3.22) the 
monosymmetric section column will buckle in pure 
Eulerian flexure, otherwise flexural torsional buckling 
failure will take place. 

If the column has two axis of symmetry then e y = 0, e y = 0 
and the characteristic buckling equation simplifies to 
(P 7 z ~N x )(Pyy -N X )(P^-N X ) = 0 (3.23) 

We find that for bisymmetric columns, the buckling 
equations become uncoupled and the three roots become 
P Z7 , Pyy and Pt}. For bisymmetric columns, the column can 
buckle in pure Eulerian flexure in the yy or zz axes, or 
pure torsional buckling mode. The smallest of the three 
roots governs the buckling mode. Critical values of each 
buckling mode is obtained in each case when m = 1. 

IV. DISCUSSIONS AND CONCLUSIONS 

The system of governing differential equations of 
buckling of an elastic column with thin walled sections 
have been solved for pinned pinned columns using the 
Fourier series method. In general, the equations were 
found to be reduced to a system of algebraic eigenvalue 
eigenvector problem. For the general use of asymmetric 
columns, the characteristic buckling equation was found 
to be identical with the Timoshenko’s solution for the 
same problem. The buckling equation for asymmetric 
columns was found as Equation (3.20). The equation is a 
cubic polynomial in N x as the unknown. The solution will 
yield three values for N x called the three critical loads as 
P C i\ > Pcr 2 anc * Per ■ AH the buckling modes in this case 

are flexural - torsional buckling modes. The critical 
buckling load, which is the lowest of the three loads 
P cr , P n - 2 and P Cf . 3 will always be smaller than the 

Euler critical flexural buckling loads Pe zz , Pryy and the 
pure torsional buckling load Pf For columns with 
asymmetric sections the buckling mode always include all 
three displacements v(x), w(x) and 8(x); hence it is a 
flexural - torsional buckling mode. 

For columns with monosymmetric cross-sections, one of 
the governing differential equations will become 
uncoupled. This results in two possible modes of buckling 
failure - pure flexural buckling failure in the yy axes (for 
columns symmetrical about the zz axes) and flexural 
torsional buckling. The characteristic buckling equation 
was found as Equation (3.22). For columns with two axes 
of symmetry, the governing differential equations of 
equilibrium become uncoupled. The characteristic 
buckling equation was found as Equation (3.23). The 


buckling equations are uncoupled and the three roots 
indicate the bisymmetric column can buckle in pure 
Eulerian flexure in the axes of symmetry or pure torsional 
buckling mode. 
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